Different mechanisms for efficient optical transmission through bilayered 
subwavelength patterned metal films 
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Light transmission through bilayered thin metal films perforated with subwavelength hole arrays 
are numerically studied based on a full-vector finite-difference time-domain approach. A variety of 
transmission peaks originating from different physical mechanisms are observed. In addition to the 
direct tunnelling and Fabry-Perot resonances, generally possessed by idealized bilayered dielectric 
slabs, the near-field localized plasmon polaritons also play important roles. They not only influence 
the direct tunnelling in a destructive or constructive way, the interactions between these localized 
plasmon polaritons on different metal films also result in additional channels which transfer optical 
energy effectively. 
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Subwavelength hole arrays in metallic thin films have 
received considerable attention recently after the experi- 
mental findings of Ebbesen in 1998 The transmission 
of light through such a nanostructured system has been 
shown to be extraordinarily efficient at specific wave- 
lengths longer than the array period 0, H, Hi- Be- 
cause of many potential applications such as near-field 
microscopy and flat-panel displays, the Ebbesen's work 
has sparked a wealth of research activities both experi- 
jliand theoretically [H E El, Q 
2ll . Moreover, various aspects of 
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the single-layer structure have been studied to exploit the 
underlying physical mechanisms, including the arrange- 
ment of the subwavelength holes [l|, Q , thickness of the 
metal film M, |T3| , polarization of the incident light f§| , 
hole shape Q and symmetry of the whole structure [2^ . 

A wealth of reasonable explanations about the Ebbe- 
sen's experiments have been presented such as the forma- 
tion of localized waveguide resonances [13] and shape res- 
onances, as well as the interference of diffracted evanes- 
cent waves [1, 0, • Among them the excitation of lo- 
calized plasmon polaritons (LPPs) is widely believed to 
playan important role 0, M, M, Q M, M, [13, M, M, 
|20| . l2l| . LPPs are the oscillations of conduction electrons 
inside the metal forced by the external electromagnetic 
waves. Because of its pure imaginary wave vector in the 
normal-to-the-interface direction, it is therefore a surface 
wave whose amplitude decays exponentially with increas- 
ing distance into each medium form the interface. This 
field confinement leads to the enhancement of the elec- 
tromagnetic field at the interface, which in turn leads 
to the extraordinary sensitive dependence of the LPPs 
properties on the surrounding dielectric conditions [23j . 
In this work, instead of studying a single-layer struc- 
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tured metal film, we investigate bilayered films with vari- 
able layer separations. We will show that at least three 
kinds of physical mechanisms contribute to the efficient 
light transmission, including the direct tunnelling, Fabry- 
Perot resonances, and additional channels formed by the 
interaction between the LPPs on different films. 

We first consider an idealized structure consisting of 
two identical dielectric slabs (with a real refractive in- 
dex 712 and thickness d) embedded into a homogeneous 
medium (with a real refractive index ni). Under normal 
incidence with wavelength A, the transmission intensity 
of an individual dielectric slab is written 



T, = 



2nin2 



1 ~i 



(1) 



with (j) = 27rn2(i/A. This equation immediately suggests 
that (f) light propagates through a homogeneous me- 
dia (assuming ni =712) without reflection; (2) the slab 
thickness d influences the optical transmission remark- 
ably. More specifically, d ~ nXjlni results in maximal 
transmission while d — {n -\- l/2)X/2n2 leads to maxi- 
mal reflection, with n being an integer. For the idealized 
double-layer structure with s being the layer separation, 
the corresponding transmission intensity is 



(2-r,)2 (T, -f)[2-2cos(a + ^)] 



T2 



r2 



(2) 



where a stands for the phase of the complex transmis- 
sion coefficient of the single slab, and = 27rnis/A. Ev- 
idently, Td — \ as long as Ts = f, in other words, an 
efficient transmission of single layer also leads to an effi- 
cient transmission of bilayered structure. We denote this 
mechanism as direct tunnelling since it only relates to 
the optical properties of single-layer slab. Eqs.Q further 
suggests Td is strongly sensitive to Lp and therefore the 
layer separation s, that is, Td oscillates between f (with 
cos{a + ip) = -I) andT2/(2-T,)2 (with cos(a-f(p) = I). 
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FIG. 1: (color online). The transmission spectrum as a func- 
tion of the layer separation s. The two dashed lines mark the 
Fabry-Perot-like resonances mainly determined by s. Ii and 
I2 indicate two additional channels induced by the interaction 
between the localized plasmon polaritons on the metal films. 
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Further assuming the phase a(A) is almost constant in 
the wavelength region of interest, the separation s con- 
sequently is proportional to the wavelength A. Because 
this mechanism is mainly determined by the separation 
s and similar to the resonant mechanisms of the Fabry- 
Perot (FP) cavities, we therefore denote it as FP reso- 
nance. It should be mentioned that the derivative ds/dX 
of different order resonance relates to the integer n and 
therefore has different value. 

The freestanding bilayered structured metal films stud- 
ied consists of two identical 500-nm-thick silver films, 
each film is further perforated by an array of square cross- 
section holes. The array period and hole side length are 
fixed to be 750 nm and 280 nm, respectively. These two 
perforated silver films are separated by vacuum with a 
variable separation s. The relative dielectric constant of 
silver is fitted by the Drude model of 



e(w) = 1.0- 



ui{lu + ij) ' 



(4) 



where cOp — 1.374 x 10 s~ is the bulk plasma fre- 
quency determined by the density of conduction elec- 
trons, 7 = 3.21 X 10^'^ s~^ represents the phenomeno- 
logical damping rate [2^. A similar single-layer silver 
film was numerically studied earlier at normal incidence, 
and a peak around 802-nm wavelength was observed in 
the transmission spectrum [l3, [111 . 

Our main result is shown in Fig.l, where the effect of 
the layer separation s on the transmission spectrum of the 



bilayered structure under normal incidence is calculated 
with the aid of a three-dimensional finite-difference time- 
domain approach We are interesting in the wave- 
length region [750, 950] nm where the far-field diffractions 
are prohibited at normal incidence because the wave- 
length is bigger than the array period. Multiple trans- 
mission peaks occur by increasing the s from 100 nm to 
980 nm with an increment of 20 nm, and three important 
properties can be observed. 

First, the subwavelength patterned metal films can be 
conceptually simplified to the idealized bilayered slabs 
studied earlier by approximating each perforated film as 
a homogeneous medium with an effective permittivity. 
Consequently physical mechanisms including direct tun- 
nelling and FP resonances are reproduced to a certain de- 
gree. More specifically, the bright column around 802-nm 
wavelength is due to the direct tunnelling, and an exam- 
ple of the electric-field distribution is shown in Fig. (2a) . 
On the other hand, the two inclined lines (marked as 
FPi and FP2 respectively) with different slopes are in- 
duced by the first-order and second-order FP resonances, 
respectively, and their typical electric-field distributions 
are plotted in Fig. (2c, 2d). Evidently, as the general char- 
acteristics of the FP resonances, different-order standing 
waves are formed between the metal films. 

Second, the strong overlap or/and repulsion between 
LPPs around two metal films lifts up the original degen- 
erate state (corresponding to s — > 00 therefore no LPPs 
interaction), and results in two additional states for ef- 
ficient light transmission [2^. Typical electric- field dis- 
tributions correspond to the higher-energy (Ii in Fig.l) 
and lower-energy (I2 in Fig.l) states are shown in Fig.(2e) 
and Fig.(2f), respectively. Evidently, the higher (lower)- 
energy state sits on the left (right)-hand side of the 
original 802-nm peak, and the corresponding wavelength 
should monotonically increase (decrease) with increasing 
the film separation s (therefore decreasing the LPPs in- 
teraction). Moreover, the LPPs interaction is so weak for 
s bigger than 900 nm that these two states almost coin- 
cide with the original 802-nm peak and hence lead to an 
enhanced transmission. On the other hand, anti-crossing 
between the second-order FP (FP2) resonance and the 
(I2) state happens around a separation of 600 nm. It 
is possibly induced by the fact that their modes greatly 
resemble each other |27| . 

Third, the separation s influences the direct tunnelling 
mechanism remarkably in the near-field region (s < 400 
nm). The transmission intensity first increases steeply 
from almost zero at s = 100 nm to roughly 80% at 
s — 280 nm (the corresponding electric-field distribution 
is plotted in Fig. (2b)), then drops to zero at s = 340 nm 
again. Although the FP effects may contribute some- 
what to this rapid variation by tuning the cos term in 
Eqs.([2]), the dominant mechanism is believed to be the 
strong overlap between the near-field LPPs around two 
metal films. The LPPs overlap constructively enhancing 
the transmission in the region s < 280 nm, while destruc- 
tively blocking the light propagation around s = 340 nm. 
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FIG. 2: (color online). The distribution of electric field ampli- 
tude corresponds to different mechanisms (see the context), 
at cross section (along the light propagation direction) with 
different separation s. (a) The direct tunnelling with s = 200 
nm; (b) The direct tunnelling with s = 280 nm; (c) The first- 
order Fabry-Perot resonance (FPi in Fig.f) with s — 350 nm; 
(d) The second-order Fabry-Perot resonance (FP2 in Fig.f) 
with s = 780 nm; (e) The higher-energy channel (Ii in Fig.f) 
with s — 660 nm; (f) The lower-energy channel (I2 in Fig.f) 
with s — 200 nm. 

Since the overlap is exponentially decreasing with the in- 
crease of the layer separation, its influence is quite weak 
for bigger s and the transmission intensity hence varies 
smoothly. 

To close the discussion, we connect our present result 
with the previous studies. It should be emphasized that 
quite similar layer-separation-transmission dependence 
was experimentally observed in the terahertz (THz) re- 
gion, and the mechanisms including both the direct tun- 
nelling and the FP resonances were found [28|. An- 
other experiment regarding a cascaded structured metal 



films sandwiched a dielectric (vacuum in our study) layer 
stressed the contributions of the FP resonance and the 
LPPs coupling [2^. Nanostructured silver multilayers 
were also experimentally investigated and the multiple 
scattering, the complicated counterpart of the FP reso- 
nance, were found to influence the optical transmission 
significantly [s^ . Few studies concentrated on the lateral 
displacement between the subwavelength hole arrays in 
different metal films. Because the LPPs overlap (inter- 
action) is sensitive to the lateral displacement, similar 
to their dependence on the vertical separation s studied 
here, the two additional states (Ii and I2) and even the 
direct tunnelling can be therefore tuned by the lateral 
displacement in the near-field region (small s) [H, [3l| . 

In summary, we studied the efficient transmission of 
electromagnetic wave through subwavelength hole ar- 
rays in bilayered metal thin films by varying the layer 
separation. Multiple transmission maximums were ob- 
served, and the underlying physical mechanisms were ex- 
ploited by comparing the subwavelength patterned struc- 
ture with an idealized bilayered dielectric slabs. The 
near-field localized plasmon polaritons at the metal films 
were found not only to modify the direct tunnelling (in 
constructive or destructive way) significantly, but also 
their interaction (overlap) form two additional channels 
for transmitting the light energy efficiently. 
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